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^^ Abstract. The interaction of trimethyl methylcyclopentadienyl platinum (MeCpPtMea) 

with a fully hydroxylated Si02 surface has been explored by means of ab initio calcu- 
lations. A large slab model cut out from the hydroxylated /3-cristobalite Si02 (HI) 

^^ surface was chosen to simulate a silica surface. Density functional theory calculations 

On were performed to evaluate the energies of MeCpPtMea adsorption to the Si02 surface. 

Our results show that the physisorption of the molecule is dependent on both (i) the 

I— I orientation of the adsorbate and (ii) the adsorption site on the substrate. The most 

O stable configuration was found with the MeCp and Mea groups of the molecule oriented 

towards the surface. Finally, we observe that van-der-Waals corrections are crucial for 
the stabilization of the molecule on the surface. We discuss the relevance of our results 
for the growth of Pt-based nanostructured materials via deposition processes such as 
electron beam induced deposition. 
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1. Introduction 

Trimethyl niethylcyclopentadienyl platinum (MeCpPtMea) is an important organometal- 
lic precursor that is widely used to deposit Pt nanostructures, not only in electron 
beam induced deposition (EBID), but also in focused ion beam deposition, chemical 
vapour deposition, laser induced chemical processing, and atomic layer deposition pro- 
cesses [IlElEliaElElIZllHllglllolIIIlIiailSlIIlIISlIIS]. 

For most of the deposition processes mentioned above, there is a lack of detailed 
understanding of the growth mechanism of nanostructures at the molecular level. While 
for instance in EBID, several Monte Carlo simulations have been performed to study 
the effects of electron energy, probe size, substrate thickness and deposit composition 
considering the primary, secondary and backscattered electrons [TTj [T8| [T9| [20t [2T| [22] , 
very little is known about the detailed geometry of adsorption configurations and the 
possible dissociation processes upon adsorption. If the various deposition processes 
are to evolve from a set of empirical recipes for producing nanostructures to viable 
scientific tools for nanofabrication, a complete understanding of the interactions between 
the substrate, adsorbed precursor molecules, and electrons in the case of EBID, is 
fundamental and must be developed jH [231123]. In Ref. [25j, a first attempt to investigate 
the adsorption of tungsten hexacarbonyl (W(C0)6) precursor molecules on a Si02 
substrate partly modified due to electron irradiation was performed by means of density 
functional theory (DFT) calculations. In the present work we consider the physisorption 
process of the more complex molecule MeCpPtMcs on a hydroxylated Si02 substrate 
within DFT with and without inclusion of long-range van der Waals (vdW) interactions. 
This study should serve as the starting point for growth of nanostructures in various 
deposition processes. Si02 is chosen as substrate in our calculations since it is one of 
the most important substrates in the semiconductor industry, and is widely used for 
instance in the EBID process. 

Physisorption properties have been studied in several systems. For example, Riniola 
et al. |26j investigated the adsorption of benzene and benzene-l,4-diol on hydrophobic 
and hydrophilic silica surfaces and found that the adsorption of the aromatic molecules 
on the hydrophobic silica surface is dictated by vdW interactions and the adsorption 
energies for the hydrophilic surfaces are almost doubled with respect to the hydrophobic 
surface due to H-bonding interactions between the substrate and the adsorbate. Mian et 
al. performed DFT calculations to investigate the catechol adsorption on silica surfaces 
and indicated that both the hydroxyls and phenylene ring of catechol contribute to its 
strong adhesion due to hydrogen bonding and vdW interactions |27j . Atodiresei et al. 
investigated the adsorption mechanism of a single pyridine molecule on Cu(llO) and 
Ag(llO) surfaces and found that vdW corrections are very important for the geometry 
and electronic structure of fiat adsorption configurations [2S| [2U] . Our results for 
MeCpPtMes on Si02 show that the adsorption of the molecule is highly sensitive to 
both (i) the molecule orientation and (ii) the adsorption site on the substrate and we 
find that only certain geometries are favored in the process. Moreover, inclusion of vdW 
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corrections are crucial for stabilizing the molecule on the surface. 

The paper is organized as follows: in the next section the methods used for the 
DFT calculations are introduced. In Section 3 the electronic properties of the precursor 
molecule, the substrate as well as the complex molecule-substrate are presented and 
discussed with special emphasis on the importance of the vdW corrections. Further, we 
discuss the relevance of our results for the deposition of nanostructures and in Section 
4 we provide a summary of our work. 

2. Methods and computational details 

Silicon dioxide (Si02) is a typical polymorphic material with structures ranging from 
amorphous to various types of crystalline phases. The /3-cristobalite crystalline phase 
and amorphous Si02 have a locally similar structure. They have similar densities 
and refractive indices, and a rather convenient cubic cell. Therefore, /3-cristobalite 
has been frequently used as a model crystal to represent the structure of amorphous 
Si02 [501 EH E21 E31 Ell ^'^^ ^i^^ ^^ ^1^° considered here as a model substrate. 

The /3-cristobalite Si02 is cubic with eight molecules in the conventional unit cell 
(Figure [I] (a)). Considering the large size of the precursor molecule MeCpPtMca, we 
need to find a compromise between computational cost and sufficiently large distances 
between adsorbates. Thus, a four-layer slab model containing 216 atoms was cut from 
bulk /3-cristobalite (space group F d2>m) with the experimental lattice constant 7.16 A 
along the [111] direction [25]. We rearranged the resulting structure by adding OH 
hydroxyl groups to the unsaturated Si atoms exposed at the surface to avoid dangling 
bond effects. The bottom O atoms resulting from the formation of the surfaces are 
also saturated with H atoms. In fact, a realistic model for a substrate prepared under 
nonvacuum conditions is hydroxylated Si02. A vacuum region of approximately 30 A 
without considering the adsorbate was placed between slab layers. The geometry is 
shown in Figure [lib) . 

Hydroxyl coverage on silica has been previously measured to be 4.9/nm^ OH groups 
at 180-200 °C and 3.6/nm^ at 300°C, independent of the structural characteristics [35] . 
Our model built by a crystalline /3-cristobalite [111] facet has a silanol density of 4.5/nm^ 
close to that observed on amorphous silica and therefore is a good representation of the 
experimental environment. 

Structure relaxation and electronic properties of the substrate and the complex 
molecule-substrate were performed in the framework of DFT by using the Vienna ah 
initio Simulation Package (VASP) [36l[37| with a projector augmented wave (PAW) basis 
set [38| |39] . Exchange and correlation effects were treated in the generalized gradient 
approximation (GGA) in the formulation of Perdew, Burke, and Ernzerhof [ID]. The 
plane wave basis set used to represent the valence electron density was cut off at 400 
eV. The Brillouin zone was sampled using a (2 x 2 x 1) F-centred mesh. Structural 
optimizations were performed by minimizing the forces on all atoms to less than 0.01 
eV/A. The two bottom layers, which are expected to be little affected by the adsorption 
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(b) 



Figure 1. (a) Bulk /3-cristobalite Si02. (b) Optimized structure of the clean surface, 
cut in [111] crystallographic direction. 
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happening at the surface, were kept fixed in the relaxation process. 

In order to compare the calculations of the adsorbed molecule-substrate system 
where periodic boundary conditions were imposed with those for the isolated molecule 
in the gas phase, we have also performed calculations for the molecule placed in a 
slightly distorted cubic cell with a cell parameter of 30 A. This lengthscale ensures a 
negligible interaction between the contents of adjacent cells and mimics the properties of 
the molecule in the gas phase. A slightly distorted cell was used to break the symmetry 
in the periodic setting. Apart from the number of k-points (a single k-point at F 
was considered for the gas phase), the same parameters were employed for (i) the gas 
phase calculations (ii) the substrate calculations and (iii) the molecule with substrate 
calculations. 

Van der Waals corrections -also called dispersion corrections- which originate 
from long-range electron interactions are not captured by the standard DFT. Several 
approaches have been developed to improve this situation |^ H3l HH |15]. In our 
calculations, we consider these effects by using an empirical DFT energy correction 
scheme (DFT-D) developed by Grinime and implemented in VASP |46] . This 
implementation has been demonstrated to be very effective for the description of 
molecular systems such as small molecular adducts, vr-stacking, and large complexes 
where dispersion corrections are relevant ^7 \ H8 | B9]. 

The dispersion corrected DFT-D energy is calculated by adding an empirical 
correction -Eoisp PS] to the Kohn-Sham energy: 

-^DFT-D = -S'KS-DFT + -S'Disp (l) 

N-l N „ij 

E-Disp = --56 2^ 2^ -^fdmp{Rij) (2) 

j=l j=i+l -"-'J 

Cq is a dispersion coefficient for the atom pair ij which is calculated by: 



CI' = \IClCi (3) 

where the coefficient C^ for atom i is given by Cg=0.05iV/pQ;*, A^ has values 2, 10, 18, 36, 
and 54 for atoms from rows one to five of the periodic table, Jp is the atomic ionization 
potential and a are static dipole polarizabilities. 

Rij denotes the interatomic distance between atoms i and j and sq is the global scaling 
factor for the functional. The term fdmp{Rij) is a damping function which is given by: 

/dmp(i?.,) = -—^QTZ^ (4) 

1 + e ^ Rr I 
The role of the damping function is to scale the force field such as to minimize 
contributions from interactions within typical bonding distances. Here, Rr is the sum 
of atomic vdW radii. Original values of all atomic and global parameters proposed by 
Grimnie are employed in the present study |46] . The summations are performed over 
all A^ atoms in the reference cell and all translations of the unit cell within a cutoff 
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Table 1. Calculated bond lengths (d) and angles ((/?) compared to the corresponding 
experimental values for MeCpPtMes. Distances are given in A. and angles are given 
in degree. The bond length dc(Cpring)-H is 1-088 A, (ic(Me)-H is about 1.099-1.103 A, 
and the bond length rfc(Cp ring)-c(Cpiing) is about 1.423-1.442 A in the calculation. 





C^Cl-Pt 


dc2--p\. 


C^CS-Pt 


'^C4-Pt 


C^CS-Pt 


Exp. |il] 
Calculated 


2.019 
2.073 


2.141 
2.075 


1.990 
2.068 


2.265 
2.376 


2.262 
2.341 




c^ce-Pt 


dc7--pt 


C^CS-Pt 


'^C4-C9 




Exp. pi] 
Calculated 


2.317 
2.321 


2.356 
2.326 


2.327 
2.346 


1.508 
1.499 






V'Cl-Pt-C2 


</3C2-Pt-C3 


V'Cl-Pt-CS 


V^C5-C4-C9 


9^08-04-09 


Expt. pi] 
Calculated 


84.1 
84.7 


86.7 
84.9 


85.9 

86.4 


127.9 
126.1 


124.0 
126.3 



radius of 30 A. The interactions over distances longer than this radius are assumed to 
give negligible contributions to -Eoisp and can be ignored. 

The relative stability of the adsorption geometries is evaluated with: 

^-E-ads = -E'system{MCPM-|-Si02) ~ -E'moleculc(MOPM, g) ~ -E'clcansurface(Si02) (5) 

where Ai^ads is the calculated adsorption energy, -E'system(M0PM+si02) corresponds to the 
total energy of the relaxed system of precursor molecule adsorbed on the substrate, 
-Emoiecuie(MOPM, g) is the total energy of the isolated precursor molecule calculated in 
the gas phase, and -Ecieansurface(SiO ) is the total energy of the clean hydroxylated Si02 
substrate. 

3. Results 

3.1. Isolated precursor molecule 



Figure pi (a) shows the optimized structure of the MeCpPtMcs molecule, which has 
one methyl substituted cyclop entadienyl ring (MeCp) and three methyl ligands (Mea) 
bound to a central Pt atom. The methylcyclopentadienyl unit is bound to the platinum 
centre atom in a pentahapto fashion with the typical platinum-carbon distances as 
found in experiment (2.262 to 2.356 A) [H]. Table hi shows selected bond lengths and 
angles obtained in our calculations compared with experimental data. The agreement 
between both sets of data is rather good. After relaxation, the distances between the 
cyclopentadienyl ring carbon and the Pt atom show changes of up to 0.11 A compared to 
the experimentally measured distances. The differences for angles between experiment 
and theory are mostly within 2°. 

In order to obtain a first view of how the charge density distributes in the molecule, 
we have calculated the electronic properties of an isolated molecule with the Gaussian 
09 code [50j by considering the standard basis sets. The structure [41j of the isolated 



Physisorption of an organometallic platinum complex on silica. An ah initio study 7 




'I 

HOMO 
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LUMO 



Figure 2. (a) Optimized structure of trimethyl methylcyclopentadienyl platinum 

(MeCpPtMea) in the gas phase with Pt in blue, carbon in gray and hydrogen in white. 
MeCp and Mea groups are marked, (b) Frontier molecular orbitals of MeCpPtMes. 



precursor molecule MeCpPtMcs was first optimized and the resulting structure is 
consistent with both experimental values and the above GGA calculation for a molecule 
in the gas phase. 

In Figure |2] (b) we show the calculated molecular orbitals for MeCpPtMea. The 
highest occupied molecular orbital (HOMO) is mostly formed by C 2p-type atomic 
orbitals from the cyclopentadienyl ring and Pt Srf-type atomic orbitals. The lowest 
unoccupied molecular orbital (LUMO) is mostly composed of a combination of 2p-type 
atomic orbitals of C bonded to Pt 5(i-type atomic orbitals. Such an orbital configuration 
indicates that a transition from HOMO to LUMO could happen by charge transfer from 
the MeCp group to the Mcs group through Pt. 
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[ SiOg surface Si02 surface SiOj surface 

(d) (e) (f) 




Figure 3. The configurations of MeCpPtMea contact to tlie surface tliat were 
considered in this work (side view), (a) MeCp ring attaching, denoted as OmcCp; 
(b) Mea group attaching, denoted as Omc-^', (c) Me from MeCp attaching, denoted 
as Omb] (d) one methyl from Mes and one methyl from MeCp attaching, denoted as 
OMe-Me; (e) two methyl from Me3 and Cp ring attaching, denoted as OMe2-Cp; (f) 
two methyl from Mcs and MeCp attaching, denoted as Omc2-McCp- 



3.2. Adsorption of MeCpPtMe^ on the hydroxylated Si02 surface 

Figure [T] shows the structure of bulk /3-cristobahte Si02 and the employed slab model in 
this study. This slab has already been relaxed, and the properties analyzed in a previous 
work |25]. With this slab, a series of adsorption geometry models based on the features 
of the MeCpPtMes molecule and of the Si02 substrate have been evaluated in order to 
screen for the most stable adsorption configuration. In Figure [3] we show the various 
configurations. If we consider the molecule as composed of Mes and MeCp ligand groups 
bonding to the central Pt, six different configurations can be considered. In three of 
them only one of the ligand groups faces the surface (configurations (a)-(c) in Figure Isl) 
and in the other three both of the ligand groups face the surface (configurations (d)-(f ) 
in Figure pi). For each configuration, a set of possible adsorption sites on the substrate 
are chosen as shown in Figures 111 (a) and (b). For sites 1-5, the orientation of the 
molecule on the surface is through either the MeCp, Me or Me^ groups (see Figure |2] 
(a) and Figure [s] (a)- (c)). For sites 6-8 both MeCp and Mes groups can face the surface 
(configurations (d)-(f) in Figure Isl). 
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(b) 



Figure 4. Classification of sites in Si02 (top view) considered suitable for adsorption 
(active sites), (a) Active sites in Si02 for one of the ligand groups (MeCp, Me or Mea 
groups, configurations (a)-(c) in Figure Isl). The sites can be classified as hollow sites 1 
and 2, bridge sites 3 and 4, and Si-OH top site 5; (b) Active sites in Si02 when both 
MeCp and Mea groups contact the surface (configuration (d)-(f) in Figure jsl, with 
hollow site 6, near bridge site 7, and top site 8. The longer border of the trapezoids 
represents the MeCp group and the shorter border represents the Mea group. 
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Table 2. Calculated adsorption energies (in eV) for various adsorption configurations 
without (denoted as -Edft) and with (denoted as _Edft-d) inclusion of vdW long 
range interaction effects. i?DFT* is the DFT energy contribution to _Edft-d- Shown 
in boldface are the energies of the most stable configurations for various molecule 
orientations. 



Orientation (Figure 3) Site (Figure 4) -Edft 



E, 



DFT-D 



E- 



O 



MeCp 



o 



Me3 



1 

2 
3 

4 
5 



1 
2 
3 

5 



DFT-D 



-0.125 

-0.128 
-0.098 
-0.178 

-0.067 



-0.465 
-0.515 

-0.452 

-0.472 
-0.283 



-0.102 
-0.109 
-0.173 

-0.068 



-0.394 
-0.438 
-0.446 
-0.496 



E- 



DFT* 



-0.138 
-0.129 
-0.073 
-0.153 
0.032 



-0.103 
-0.128 
-0.047 
-0.030 



Omc 


1 


-0.077 


-0.327 


-0.063 


2 


-0.111 


-0.359 


-0.013 


Omc-Mc 


7 


-0.125 


-0.526 


0.148 




6 


-0.101 


-0.477 


-0.108 


Olvlea-Cp 


7 


-0.148 


-0.669 


-0.149 




8 


-0.184 


-0.596 


-0.098 




6 


-0.111 


-0.503 


-0.129 


OlvIea-MeCp 


7 


-0.154 


-0.495 


-0.089 




8 


-0.078 


-0.650 


0.000 



Table [2] lists the calculated adsorption energies where the structure relaxation and 
the energy calculations were performed without (denoted as -Edft) and with inclusion of 
vdW long range interaction effects (denoted as Edyt-h)- We have also included a third 
column in the table denoted by -Edft* which contains the DFT energy contributions to 
-Edft-d (i-e. vdW energy corrections have been substracted from i^DFT-D)- 

In order to analyze these energies, we have plotted in Figure |5] the binding energies 
of MeCpPtMcs on the Si02 surface as a function of the active adsorption site without 
(DFT) and with inclusion of vdW corrections (DFT-D) for all considered configurations. 
Also shown is the decomposed DFT contribution in DFT-D (DFT*). Within the DFT 
calculation, the adsorption site and the orientation of the molecule are found to be 
equally important for the adsorption energy. With inclusion of vdW corrections, the 
adsorption energy is primarily determined by the orientation of the molecule and, to a 
lesser degree, by the chosen adsorption position of the molecule on the substrate. For 
different molecular orientations, the DFT-D adsorption energy varies by about 310 meV; 
meanwhile, for a given orientation, changes in the adsorption position on the substrate 
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Figure 5. (a) Calculated DFT, DFT-D and DFT* adsorption energies; lines are a 
guide to the eyes, (b) Distances between Pt atom and the topmost Si plane. 



lead to variations in the adsorption energy between 32 and 232 meV. The most stable 
configuration is significantly modified when vdW corrections are included. In fact, the 
effect of vdW forces to the weak bonding of MeCpPtMes to a fully hydroxylated Si02 
surface is to bring the molecule closer to the surface. 

If we consider the centered Si from the substrate as the topmost Si plane location, 
the distance between the Pt atom and this Si plane for different configurations in the 
DFT-D calculation is shown in Figure [S] (b). As we can see, the difference between 
various configurations follows a similar trend as the adsorption energy calculated in 
DFT-D (Figure p^ (a)). The vdW correction brings the adsorbate molecule closer to 
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(b) 



Figure 6. The most stable structure found of MeCpPtMea adsorbed on Si02 surface, 
with (a) for side view and (b) for top view. Distances in (a) are given in unit of A, 
relevant surface silanols are numbered in blue. 
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the substrate when both MeCp and Mes groups face the surface, thus enhancing the 
interaction at the interface. According to Xue et al. [41j, it is difficult to see how the 
platinum can contact the surface in the ground state of the molecule so that some 
dislocation of the ?75-MeCp ring to lower hapticity is required before the platinum atom 
of the hydrocarbon complex may interact with the surface atoms. In our case, the 
distance between Pt and the substrate in the geometry where both Mea and MeCp 
groups face the surface is lower than for other orientations (Figure [s]), indicating that 
this configuration enhances Pt-substrate interaction strength. 

As shown in Table |2] and Figure |5| the molecule prefers to stay on the surface with 
both Me3 and MeCp ligand groups oriented towards the surface (see Figure IS] (d)-(f)) 
while the orientation with only one of the ligand groups oriented towards the surface 
is about 154 to 310 meV less stable (see Figure [3] (a)-(c)). Note, that the 0Me2-MeCp 
configuration on site 8 is only 19 meV higher in energy, making this structure a close 
runner-up to the optimal configuration 0Me2-Cp on site 7; in fact both structures are 
rather similar and the energies are difficult to distinguish with the precision that we 
could achieve in our calculations. 

In Figure |6] we show the most stable adsorption configuration calculated with 
inclusion of vdW corrections found in our work with the MeCp group oriented 
perpendicular to the surface and coupled with one edge of the aromatic ring to the 
substrate (Omc2-Cp)- After relaxation of the molecule-substrate complex, we observe no 
significant structural deformation of the precursor molecule. The optimized distances 
between the topmost O atoms on the substrate and the closest atoms of the MeCpPtMcs 
molecule are shown in Figure |6] (a). The spacing is significantly larger than typical 
covalent or hydrogen bond distances. The structural properties indicate that the 
precursor molecule is indeed physically adsorbed on the fully hydroxylated Si02 surface 
without formation of covalent bonds with the substrate. 

From the calculated binding energies at different sites, we can estimate the surface 
diffusion barriers, which determine the diffusion ability along different directions [51j. 
At low coverage, the precursor molecule binds to the surface with both Me3 and MeCp 
ligand groups. The orientation 0Me2-Cp prefers site 7 while 0Me2-MeCp prefers site 8, 
with a large thermodynamic diffusion barrier of 73-155 nieV along the [211] direction. 
At high coverage, the molecule can also bind to the surface with only one of its ligand 
groups, in the orientations OmcCp and Oiviea, with small diffusion barriers. 

After adsorption, the C-C and C-H bonds of the molecule are slightly weakened, 
whereas the Pt-C bonds of the Pt-Me3 group are slightly shortened. The distance 
between Pt and the MeCp ring are significantly decreased with the length between Pt 
and C atoms of the MeCp ring shortened by 0.068-0.133 A. The surface region far 
from the adsorbate shows almost no changes. However, the 0-H bond of the surface 
silanol and some Si-0 bonds close to the adsorbate are weakened after adsorption. This 
observation indicates that the precursor molecule is loosely bonded to the substrate. 
We found that Pt dyz and ^^2 orbitals and the C Py orbitals (the orbital notation is 
given in the local reference frame of the molecule where the xy plane is defined parallel 
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to the MeCp group) play an important role in the interaction between the MeCpPtMea 
adsorbate and the surface. 

The orbital wave function of the HOMO precursor molecule (Figure |2](b)) indicates 
that the 0Me2-Cp orientation has the largest charge density near the substrate orbitals 
and can thus lead to more interaction than other orientations. This explains why this 
is the most favorable configuration on the surface. In contrast, the small contribution 
of the Me from the MeCp ligand to the HOMO orbital explains why the configuration 
where this ligand group faces the substrate corresponds to the most unstable adsorption 
configuration (the most unfavorable adsorption energy found in this work (see Table [2]). 
In previous studies on adsorption of molecules containing conjugated vr electrons in 
aromatic rings like cyclopentene (c-CsHg) on Ni (111) |52], cyclop entadienyl anion (c- 
CsHs") on Ni (111) [53] or pyridine on Cu (110) and Ag (110) surfaces [28], the geometric 
orientation in equilibrium was found to be the one where the cyclic hydrocarbons 
orient nearly parallel to the surface. In those cases the vr systems contribute to the 
interaction between adsorbate and substrate. This is not the situation in the present 
more complicated MeCpPtMea precursor molecule. In our case, the pz states of C from 
the Cp ring are located at about -5 eV deep down in the valence band. 

In Fig. u^we show the total density of states (DOS) of the preferential configuration 
of MeCpPtMea on Si02 as well as the partial contributions of the adsorbate. Results 
are shown for calculations including vdW corrections. The highest molecular level of 
the MeCpPtMes adsorbate lies in the gap of the substrate, which pins the Fermi level 
close to the valence band. The largest contribution at the Fermi level are Pt and MeCp 
states. We observe that the DOS of the Si02 substrate shows only minor modifications 
with respect to the DOS of the free Si02 substrate (compare with Fig. 6 (a) from Ref. 
|25|). These results confirm the weak interaction between the adsorbate molecule and 
the Si02 substrate. 

We can now relate our results to some deposition processes like EBID. Electron 
irradiation of adsorbed MeCpPtMes produces hydrogen and methane as the only volatile 
products [3], and the methane could be associated with cleavage of Pt-CHs (methyl 
from Mes) or C-CH3 (methyl from MeCp). Figure u\ shows that hydrocarbons that 
contribute to the valence band states near the Fermi level are mainly located in the 
Mes group and the Cp ring of the molecule. This indicates that the methyl on the 
Mes group is more active than that of the MeCp group, therefore theoretically they 
will be easier to dissociate when external energy is applied. This is in agreement 
with recent experimental work where the same precursor molecule decomposition on 
a gold substrate in an EBID process has been investigated with X-ray photoelectron 
spectroscopy (XPS) and mass spectrometry (MS) |2]. The authors of Ref. [3] find 
that trimethyl-methylcyclopentadienyl platinum adsorbed on Au substrates at 180 K 
undergoes electron stimulated decomposition mediated by Pt-CHs bond cleavage instead 

of C-CH3. 
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Figure 7. DOS of the preferential configuration of MeCpPtMes on Si02. Shown 
are (a) the total DOS and (b) the partial contributions of the adsorbate. The Fermi 
energy is set to zero. 



4. Summary 

Our calculations based on DFT show that the MeCpPtMes molecule is weakly bonded 
to the fully hydroxylated Si02 surface without deformation. The calculated adsorption 
energy is found to be -0.669 eV where the molecule is located at a nearly hollow site, 
and the distance of Pt to the topmost -OH surface silanol group is 2.87 A. Contrary 
to the adsorption of molecules containing conjugated tt electrons in aromatic rings, the 
preferential configuration of MeCpPtMcs corresponds to having the MeCp and Mes 
groups oriented towards the surface. We find that the adsorption is more dependent 
on the orientation of the adsorbate with respect to the surface than the adsorption 
site. We also observe that the optimal configuration of the molecule corresponds to a 
minimum in the Pt-substrate distance, indicating that the interaction between platinum 
and substrate is maximized in the lowest energy configuration. Moreover, inclusion of 
vdW corrections contribute to the stabilization of the molecule on the surface. Our 
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calculations provide a theoretical insight on the interfaces between MeCpPtMcs and a 
Si02 substrate, which is essential for the growth processes of deposits. 
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